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• Background and Aims The root apical meristem (RAM) is the plant stem cell niche which provides for the 
formation and continuous development of the root. Auxin is the main regulator of RAM functioning, and 
auxin maxima coincide with the sites of RAM initiation and maintenance. Auxin gradients are formed due to 
local auxin biosynthesis and polar auxin transport. The PIN family of auxin transporters plays a critical role 
in polar auxin transport, and two mechanisms of auxin maximum formation in the RAM based on PIN-mediated 
auxin transport have been proposed to date: the reverse fountain and the reflected flow mechanisms. 

• Methods The two mechanisms are combined here in in silico studies of auxin distribution in intact roots and 
roots cut into two pieces in the proximal meristem region. In parallel, corresponding experiments were performed 
in vivo using T)K5\\G¥V Arabidopsis plants. 

• Key Results The reverse fountain and the reflected flow mechanism naturally cooperate for RAM patterning and 
maintenance in intact root. Regeneration of the RAM in decapitated roots is provided by the reflected flow mech- 
anism. In the excised root tips local auxin biosynthesis either alone or in cooperation with the reverse fountain 
enables RAM maintenance. 

• Conclusions The efficiency of a dual-mechanism model in guiding biological experiments on RAM regener- 
ation and maintenance is demonstrated. The model also allows estimation of the concentrations of auxin and 
PINs in root cells during development and under various treatments. The dual-mechanism model proposed 
here can be a powerful tool for the study of several different aspects of auxin function in root. 

Key words: Auxin response, root apical meristem, patterning, reverse fountain, reflected flow, mathematical 
model, Arabidopsis thaliana. 



INTRODUCTION 

In higher plants post-embryonic growth is maintained by the 
activities of meristems, which comprise stem cells and their 
derivatives, serving as the source of all differentiated cell 
types. The root apical meristem (RAM) is established at the 
earliest stages of embryogenesis and is localized to the root 
apex after germination (Dolan et al., 1993). Accumulating evi- 
dence suggests that the plant hormone auxin has the primary 
role in RAM establishment and patterning (reviewed by Jiang 
and Feldman, 2005). The distribution of auxin in the root 
apex as inferred from the activity of auxin-responsive promoter 
elements DR5 and DII has a concentration maximum in the 
RAM, namely in the quiescent centre (QC) and root cap initials, 
with a decreased level in the root cap (Sabatini et al, 1999; 
Brunoud et al, 2012). The existence of an auxin maximum in 
the QC was also shown by combining fluorescence-activated 
cell sorting and highly sensitive mass spectrometry methods 
(Petersson et al, 2009). This maximum of auxin concentration 
forms at the earliest stages of RAM initiation, persists at the root 
tip in development and is restored from provascular cells after 



RAM structure damage (before RAM regeneration) (Sabatini 
etal, 1999; Xu a/., 2006). 

Auxin gradients have been shown to be formed due to local 
auxin biosynthesis and/or directional cell-to-cell polar auxin 
transport (Sabatini et al., 1999; Benkova et al, 2003; Friml 
et al, 2002; Stepanova et al, 2008; Ikeda et al, 2009; 
Petersson et al, 2009; Kim et al, 2011). Auxin transport is 
driven by the asymmetric distribution of different auxin mem- 
brane carriers, including AUX/LAX protein family members, 
p-glycoprotein ABC transporters and PIN family efflux car- 
riers (Bennett et al, 1996; Friml et al, 2003; Blilou et al, 
2005; Geisler and Murphy, 2006). PIN family efflux transpor- 
ters were shown to have a rate-limiting role in auxin efflux and 
direct auxin flow in plants (Petrasek et al, 2006; Wisniewska 
et al, 2006). Compared with ABC transporters, the action of 
PINs is auxin-specific and more sensitive to polar auxin trans- 
port inhibitors (Petrasek et al., 2006). Pin mutants resemble 
phenotypes of the wild-type plants treated by auxin transport 
inhibitors in more aspects than do mutants of ABC transporters 
(reviewed by Tanaka et al, 2006). In Arabidopsis, the 
PIN-FORMED (PIN) proteins are encoded by functionally 
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redundant genes so the defects in pin mutants are masked by 
ectopic activity of the remaining PIN genes (Bhlou et al, 
2005; Vieten et al, 2005). In root development, phenotypic 
defects were revealed only for triple pinl pin3 pin4 mutants 
having very short root and for quadruple pinl pin3 pin4 pin7 
mutants, which, depending on the ecotype, either are embryo 
lethal with pronounced root pole abnormalities or develop 
into seedlings with no root or a non-functional root (Friml 
et al, 2003). On the other hand, loss of ABC or AUX/LAX 
transporters does not appear to strongly impact early develop- 
mental processes (reviewed by Ugartechea-Chirino et al, 
2010; Zazimalova et al, 2010). This provided additional 
evidence for a critical contribution of PIN proteins to the 
auxin distribution pattern in the root. In particular, the PINl, 
PIN3, PIN4 and PIN7 proteins provide for a continuous acrop- 
etal (rootward) auxin flow via the vascular system to the QC 
cells (Friml et al, 2002; Blilou et al, 2005). PIN3 and PIN7 
are also involved in the lateral redistribution of auxin in the 
root cap. PIN2 proteins mediate basipetal (shootward) auxin 
transport from the root tip via the epidermis as well as rootward 
auxin transport in root cortex. Here we use the rootward/shoot- 
ward directional nomenclature defined by Baskin et al. (2010). 

The assumption that PIN family efflux carriers are the main 
contributors to the formation of the auxin distribution pattern 
in the root was central to two different hypotheses regarding 
the mechanism responsible for auxin distribution within 
the root: a 'reverse fountain' (Swarup and Benett, 2003; 
Grieneisen et al, 2007) and a 'reflected flow' (Likhoshvai 
et al, 2007; Mironova et al, 2010). According to the 
concept of the reverse fountain, the rootward and oppositely 
directed shootward auxin flows are coordinated by lateral 
auxin redistribution processes to generate and maintain an 
auxin distribution in the root tip. This mechanism, whereby 
tissue patterning predetermines the morphogen distribution, 
could be characterized as a 'structural mechanism' due to 
the pre-existing structure of PIN-localized expression levels. 
The 'reflected flow' hypothesis could be referred to a 'pos- 
itional information' mechanism, where, by contrast, the 
morphogen distribution determines tissue patterning and ana- 
tomical structure (Mironova et al, 2010; Wolpert, 2011). 
The latter mechanism is based on the auxin-dependent regula- 
tion of auxin flow: low auxin concentrations activate the tran- 
scription of PINl genes, whereas high concentrations induce 
degradation of PINl proteins. The two mechanisms were 
tested in two-dimensional (2-D) mathematical models which 
describe auxin redistribution by diffusion and active transport 
in a cell layout located at a longitudinal cut of root. 

We suggest that the reverse fountain and the reflected flow 
mechanisms are complementary in root development. In par- 
ticular, only the reflected flow mechanism operates at the 
very early stages of root development. At later developmental 
stages, an anatomical structure forms and provides for the 
functioning of the reverse fountain mechanism, which allows 
for more robust maintenance of the auxin maximum in the 
RAM. However, the reflected flow mechanism does not dis- 
appear, revealing itself if RAM structure is disrupted. To test 
this hypothesis on the complementarity of the two mechan- 
isms, we combined them in a dual-mechanism 2-D mathemat- 
ical model. In numerical experiments it was shown that the 
reverse fountain and the reflected flow mechanisms could 



naturally cooperate due to auxin-regulated expression of 
auxin transporters. We also investigated the distribution of 
auxin response in DR5-GFP seedlings with decapitated roots 
and root tips excised from these seedlings. In decapitated 
roots the reverse fountain is destroyed and in silico experi- 
ments showed that the reflected flow mechanism alone may 
provide for RAM restoration after the root tip is cut off. 
However, experimentally all excised roots, including the smal- 
lest roots with the least tissue left above the QC, were able to 
maintain the positioning of the auxin maximum for more than 
4 d. This could not be explained by the functioning of either 
the reverse fountain or the reflected flow mechanism, alone or 
in combination. We tested the role of auxin synthesis in mainten- 
ance of an auxin maximum. In these numerical experiments we 
showed that ubiquitous auxin production at a small but non-zero 
rate provides for an auxin distribution pattern that agrees with 
the experimentally recorded auxin response dynamics in the 
excised root tips. Thus, both computer simulations and experi- 
mental validation indicate that auxin biosynthesis is a third 
force of equal importance to the reverse fountain and reflected 
flow mechanisms for auxin patterning in the RAM. 

MATERIAL AND METHODS 

Plant material and root excisions 

DR5::GFP seeds were kindly provided by Alexis Peaucelle 
(INRA Centre de Versailles-Grignon, Versailles, France). 
Arahidopsis seedlings were grown under a 16 : 8-h light- 
dark cycle at 18 °C on vertically oriented Petri dishes with 
half-strength MS agar medium containing 1 % sucrose. 

The root cut experiments were performed as shown in Fig. 3 
below. The root tips of seedlings at 3 d after germination (dag) 
were dissected manually by a steel blade under a stereomicro- 
scope. An incision was made to the proximal meristem, the 
exact position varying from plant to plant. Over the next 4 d (6, 
24, 48, 72 and 96 h) after dissection changes in green fluorescent 
protein (GFP) expression were fixed in both the decapitated root 
and the excised root tip using an LSM 510 META laser scanning 
microscope or a Zeiss Axiolmager Ml. Images were processed 
using ImageJ and AxioVision software. RAM regeneration was 
confirmed by tip morphology and its gravitropic response. 

2-D dual-mechanism mathematical model 

Model assumptions. In the present model we consider processes 
occurring in the RAM. To simplify the model, only three auxin 
transporters are considered: PINl, PIN2, and PIN3. PINl 
protein mediates acropetal (rootward) auxin flow in the central 
cylinder cells; PIN2 is expressed in the epidermis and regulates 
basipetal (shootward) and lateral inward transport; and PIN3 
regulates lateral outward auxin redistribution in the root tip 
where auxin accumulates. We ignore the redundant mechanisms 
for PIN2 and PIN3 to mediate rootward auxin flow and for PINl 
to mediate lateral inward auxin transport. 

Discussion of PIN protein polarization mechanisms are 
beyond the scope of the present work. It is assumed that all 
the PIN protein synthesized in a cell is localized to the cell 
wall. The localization of PIN proteins is different for the two 
cell types. PINl is located only at the basal (rootward) side 
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of the central cylinder cell, PIN2 proteins at the apical (shoot- 
ward) and lateral (inward) cell sides of the epidermis, and 
PIN3 is expressed and localized ubiquitously (Fig. IB). 

The concentration of PIN proteins in a cell depends on auxin 
concentration in a cell and is described according to the experi- 
mental data of Vieten et al. (2005). Auxin activates PINl-4 
and PIN7 transcription at relatively low concentration, 
whereas at high concentrations, auxin mediates increased deg- 
radation of PINl, PIN2 and PIN7 proteins. Data regarding the 
effective concentrations of auxin for PINl, PIN2 and PIN3 
regulation from Vieten et al (2005) (Fig. IC, D) were used 
for model parameter estimation (see below). 

The DR5 reporter was demonstrated to reflect free auxin 
concentrations between 10~^ and 10""^ m (Ulmasov et al, 
1991 \ Sabatini et al, 1999). In the model analysis we 
assume that the auxin concentration in a cell correlates linearly 



with the auxin response observed by DR5::GFP. This allows us 
to compare in vivo and in silico results. 

Model description. In the model the dynamics of auxin concen- 
tration [a] distribution and expression of its transporters 
[PINl], {PIN2] and {PIN 3] in the root meristem are described 
based on the interactions of a limited set of intracellular 
processes, as follows. 

Subsystem 1: auxin flow to the root meristem 
The rate of auxin flow (Va) into the meristem from the root 
elongation zone (shoot-derived or rootward auxin flow) is 
described by: 

Va = a (1) 
where a is the intensity of auxin flow. 



A Epidermal 

Central cylinder ^ 

1 i- 



M 



it 



N 



B 



PIN1 



PIN2 



■ ■ ■ ■ 




r r r r > 


■ 






■ 










r r 


■ 




r / r 



PIN1 



PIN2 



PINS 




Low 
auxin 
> 





auxin 

--->u 



High 
auxin 



Low 
auxin 
> 

High 



Low 
auxin 



->J2f 
->0 
->0 



D 




0-01 jlM 1 |LIM 10 jlM 

Auxin concentration, [a] 



Fig. 1 . Schematic representation of the main model assumptions. (A) The cellular layout of the 2-D dual-mechanism model, which represents a longitudinal cut 
of the root meristem. Only two cell types were specified in the M x rectangular layout: epidermal tissue (grey; j = 1,2 and j = M- 1, M) and central cylinder 
(green; j = 3, . . . , M-2). Layer / = 1 is located at the root end; i = N is at the border of the meristematic zone. The auxin transport routes by PINl and PIN2 
proteins are presented by red and blue arrows, respectively. For simplicity, PINS, which is considered to transport auxin without specific direction, is not shown. 
(B) PIN localization within the cells of different types. PINl is located basally in a central cylinder cell; PIN2 has basal-lateral expression in the 'epidermal' 
tissue; and PINS may be ubiquitously expressed along the cell membranes of both tissues. (C) The mechanisms of auxin-regulated PIN expression considered in 
the model according with the data from Vieten et al. (2005). (D) Qualitative plot of PIN protein expression versus auxin concentrations based on the mechanisms 

(C) and experimental data (Vieten et al, 2005). 
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Subsystem 2: auxin biosynthesis The rate of degradation for PIN proteins Vd,piN is defined by 

We examined two variants of auxin synthesis (with the rate the rate equation: 
Vs a) in a root: no synthesis (2a) or auxin is synthesized in a 

ceil at a low rate (2b). V^^pin = K^,pin • [PIN] 'f{[a]), (7) 



Vs,a = 0 (2a) 

Vs,a = ^s,a (2b) 

where K^ ^i is the constant of auxin synthesis. 
Subsystem 3: auxin degradation 

The rate of irreversible loss (degradation) of auxin, y^^a, 
through its utilization or migration from the modelled region 
is defined as 

Vd,a = ^d,a • W (3) 

where K^^^^ is the rate constant for auxin degradation and [a] is 
the concentration of auxin. 
Subsystem 4: auxin diffusion 

We ignore auxin diffusion within the apoplast; instead, 
auxin is modelled as moving directly from a cell to its neigh- 
bours by means of passive diffusion and active transport. The 
rate of auxin diffusion from one cell to another is described as: 

Vdiff,a = • [a] (4) 

where D is the diffusion rate constant. 

Subsystem 5: auxin active transport mediated by PIN 
transporters 

The rate of active transport from one cell to another is 
described by the following mass action equation: 

Vu,a = Ko,PiN • [a] . [PIN] (5) 

where [PIN] denotes concentration of the corresponding trans- 
porter (PINl, PIN2, PIN3) and Kq^pin are the constants of 
active transport rate via the corresponding PIN (/^Tq^pini for 
PINl and Kq. ,piN3 for PIN3). As PIN2 mediates both shootward 
and inward auxin flows in epidermis, ytr,a is divided between 
the basipetal and lateral directions in the ratio 

^0,bas,PIN2 • Ko,lat,PIN2 = 7:3. 



where K^^pin are the degradation constants of the correspond- 
ing PIN proteins (PINl, PIN2, PIN3). As auxin regulates PINl 
and PIN2 degradation, the non-linear part of the degradation 
rate/([(2]) was defined as: 

/ [a] 

/(M) = l + 5™^ (8) 

where q3,^piN are the thresholds of auxin- dependent PIN deg- 
radation, and hpiN are the coefficients that define non-linearity 
of auxin-regulated PIN inhibition. The 8pi^ coefficients define 
if the PIN degradation is regulated by auxin (Spj^j^ — ^pin2 — 
1, ^IN3 = 0)- 

Model construction 

Elementary subsystems 1-7 were used as the blocks for the 
2-D dual-mechanism model construction according to the 
model configuration (Fig. 1). The rectangular cell layout con- 
sists of M layers ( j = 1, . . . , M) of A/^ ceUs each (/ = 1 . . . N), 
where rows j = 3, . . . , M - 2 correspond to the central cylinder 
cells and rows j= 1, 2, M- 1, M correspond to epidermal 
tissues. Cells in layer / = N are the most proximal cells of 
the root meristem, whereas cells in layer / = 1 are the root 
end. For all cells in the layout subsystems 2 and 3 are consid- 
ered (in the basal variant of the model there is no auxin syn- 
thesis, eqn 2a). Also in all cells subsystems 4 and 5 operate 
taking into account the neighbours to which auxin moves by 
means of diffusion and active transport. In a cell of the 
central cylinder, subsystems 5-7 are integrated for PINl and 
PIN3, but for an epidermal cell subsystems 5-7 are defined 
for PIN2 and PIN3. Subsystem 1 is defined only for the / = 
N cells of the central cylinder. The system of ordinary differ- 
ential equations for the 2-D dual-mechanism model is 
assembled by summation of all rates (1-7) over the cellular 
layout (see Supplementary Data for details). 



Subsystem 6: synthesis of PIN proteins 

The synthesis rate of PIN proteins Vs,pin was approximated 
by the following Hill function taking into account the positive 
regulation by auxin: 



' s,PIN 



[a] 

<ll,PIN 



s.PIN ' 



1 + 



a]_y 



(6) 



Q2,PIN/ 



where K^^pm are the synthesis rate constants for the corre- 
sponding PINs (PINl, PIN2, PIN3), qi^piN are the thresholds 
of auxin-dependent activation of PIN synthesis, q2,piN are the 
thresholds for saturation of auxin-dependent PINl synthesis 
and Spijsi are the Hill coefficients. 
Subsystem 7: degradation of PINs 



Parameter estimation 

We used the experimental data of Vieten et al (2005) as a 
basis for parameter estimation. The effective auxin concentra- 
tions for PIN upregulation revealed in the experiment are 
presented in Supplementary Data Table SI. We estimated 
the parameters for Hill functions so that the maxima of PIN 
expression in silico coincide with the effective auxin concen- 
trations measured by Vieten et al. (2005). The other para- 
meters were adjusted so that the patterns of PIN expression 
and auxin distribution qualitatively match the experimental 
data (Supplementary Data Table S2). 



Numerical simulations 

The non-linear system of equations of the model was 
integrated using the MGSModeller software package 
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(Likhoshvai et ah, 2001) on high-performance supercomputer 
cluster SSCC (Siberian Supercomputer Center). The dual- 
mechanism model calculation and analysis were repeated in 
the Matlab system (see Supplementary Data). Depending on 
the model configuration and the set of parameters, the calcula- 
tions were done until a stationary solution was achieved 
(on average approx. lO"^ iterations). 

The model simulation results were visualized using the 
Paraview program (Cedilnik et al, 2006). 

RESULTS 

Cooperation of the reverse fountain and reflected flow 
mechanisms in the RAM 

To investigate cooperation of the reverse fountain (Grieneisen 
et al, 2007) and the reflected flow (Mironova et al, 2010) 
mechanisms in the root meristem the 2-D dual-mechanism 
model was created as follows (see also Material and 
methods, Fig. 1). The dual-mechanism model is based on 
the 2-D reflected flow model structure (Mironova et al, 
2010), to which the processes of PIN2 and PIN3 expression 
were added to provide for shootward and lateral (outward/ 
inward) auxin flows in the root meristem. The dual-mechanism 
model differs from the reverse fountain model in considering 
the dependence of PIN (PINl, PIN2 and PIN3) expression 
on auxin concentration in a cell. Semi-quantitative estimation 



of the dependences were made by using the data from 
Vieten et al, (2005) (Supplementary Data Table SI), who 
showed that the amounts of PIN transcripts and PIN proteins 
in the root meristem are strongly correlated with the concentra- 
tion of exogenous auxin. 

We adjusted the parameter values for the dual-mechanism 
model such that patterns of PIN expression and their mode 
of auxin dose dependence fit the experimental data (Fig. IC, 
D; Supplementary Data). The dual-mechanism model with 
this set of parameters and initial uniform distribution of 
auxin and PIN protein concentrations in the cell layout 
yields a stationary solution with auxin and PIN protein patterns 
matching the experimentally observed ones (Fig. 2; 
Supplementary Data Video SI). 

The values and the dimensionality of the effective auxin con- 
centrations for upregulation of PINs revealed by Vieten et al 
(2005) were used to estimate the values of variables in the dual- 
mechanism model. Auxin concentration was estimated as 04- 

1 |JLM in a vascular cell, 10-11 (jlm in the QC and 8-10 |jlm in 
a columella cell; PINl in a vascular cell was estimated as 

2 |JLM, PIN2 in a epidermal ceU as 04-1-2 |jlm and PIN3 in a 
columella cell as 04 -0-65 |jlm. Note that the estimations may 
vary significantly by taking into account other processes that 
affect auxin transport. For example, the model does not repro- 
duce high PIN3 expression in the central cylinder (Fig. 2C), 
where PIN3 may be upregulated by other mechanisms. It has 
been shown that cytokinin- auxin crosstalk antagonizes BRX 



A PIN1 B PIN2 




C PINS D Auxin 




Fig. 2. The 2-D dual-mechanism model stationary solution. For each variable, a schematic representation of the expression pattern is given compared with the 
steady-state model solution. (A) PINl expression; (B) PIN2 expression; (C) PINS expression; (D) DR5 expression revealed in vivo and auxin distribution calcu- 
lated in silico. 
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Fig. 3 . In vivo and in silico experiment on RAM regeneration after root cutting. (A) Scheme of the in vivo experiment on root cutting. (B) Simulation scheme of 
root cut in the model. (C) Regeneration of the auxin response maximum in the decapitated root. (D) Correspondence between the in vivo and in silico results on 
changes in auxin response measured along the central root axis during regeneration of the decapitated root. Auxin (response) distribution is shown 2 d after 
cutting. Left: the dynamics of the DR5 signal along the central root axis of the regenerating RAM. Right: auxin distribution in the central cylinder layer j = 4 
obtained at different time points of the model solution as the dual-mechanism model approached the stationary state. (E) Scheme of the model results on simulation 
of root regeneration after decapitation. Three different scenarios of the new root tip development were described dependent on the size of the PM that was left in the 

decapitated root both in vivo and in silico. Scale bar =100 \xm. 



activity in the stele, which may affect PINS expression there 
(Scacchi et al, 2010). 

The mechanism of auxin distribution self-organization 
found in the resulting stationary solutions is as follows. At 
the first early phase, an auxin maximum is generated in the 
central cylinder cell array at a distance from the root end 
under the reflected flow mechanism (Mironova et al, 2010). 
As a result, a zone of high auxin in the root tip is organized, 
which determines the peak of PINS synthesis. At the second 
phase, PINS-mediated auxin redistribution is switched on in 
the root tip and auxin moves laterally from rootward flow 
and enters the epidermis where it is further involved in 
PIN2-based shootward flow. As PIN2 is localized also on 
the lateral internal cell sides in the epidermis, it redistributes 
auxin back to the rootward flow, finally completing the 
structure of the reverse fountain at the third phase. 

Thus, our simulations demonstrate that the reflected flow 
and the reverse fountain may naturally cooperate in providing 
for establishment and further maintaining of the auxin pattern 
in the root tip. 

Auxin redistribution after root cut 

To estimate the advantage given by the cooperation of the 
two mechanisms, we made in vivo and in silico experiments 



on root cut. Root tip isolation abolishes the reflected flow 
mechanism in this tip, whereas cutting off the root tip 
damages the reverse fountain structure in the decapitated root. 

In vivo experiments on root cuttings. We cut roots of DR5::GFP 
seedlings (S dag) in such a way that both cuttings have part of 
the proximal meristem (PM) (Fig. SA). Over the next 4 d we 
monitored reporter activity in both the decapitated root and 
the excised root tip. 

In the tip of the decapitated root we observed a restoration of 
the auxin response peak followed by meristem regeneration. 
If the decapitated root contained less than half of the initial 
PM, the new root end shows high DR5 activity over the first 
2 d, but the RAM never regenerates from it. If it contains 
between one-half and approx. four-fifths of the initial PM, 
the cut edge forms callus-like outgrowth first, followed by re- 
generation of the root meristem from the outgrowth periphery. 
The outgrowth cells have moderate and non-uniform DR5 
activity, this activity increases in the regenerating meristem. 
Finally, if the new root end contains more than four-fifths of 
the PM, but still does not include the QC, root meristem regen- 
eration occurs after S d without the callus stage. In this case, 
the dynamics of DR5 activity were similar to those observed 
by Xu et al. (2006) after QC laser ablation. During the first 
12 h DR5 activity increased at the cut edge (Fig. SC). The 
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region of DR5 activity expands in time so that the new auxin 
response peak was formed proximally to the cut edge in first 
24 h. Auxin response maximum positioning and subsequent 
RAM regeneration occur over the next 48 h. 

In the excised root tip the auxin pattern was maintained and 
the RAM remained viable for a long period of time. Although 
DR5 activity decreased over time (Fig. 4A-E), the signal was 
observed in the root tips at least for 4 d after cutting. The auxin 
response persists even in the shortest isolated tips without any 
vasculature. 

In silico experiments on RAM regeneration. We simulated the 
auxin maximum restoration processes after root tip cut in the 
following manner (Fig. 3B). We first calculated the stationary 
auxin concentration [a]i j in the dual-mechanism model at M = 
8 and N =?>5 and the set of parameters (Supplementary Data 
Table S2). The auxin distribution in the solution has a 
maximum in the 4th cells of the central cylinder (Fig. 2D). 
We than simulated the cut off root tip with a length of k 
cells and modelled the cell layout at A/^' = 35 - k. The calcula- 
tions were made with the same set of parameters 
(Supplementary Data Table S2) but the initial data were 
[a\j = [a]i+kj and [PIN% = [PIN]i+k,p where / = 1 . . . 35 - /: 
and the variable [PIN] was [P/W], {PIN2^ or \PIN3\ The 
configuration of the model with the specified auxin and PIN 
distribution simulates the decapitated root tip. Note that in 
the initial model the length of the PM was N-A = 31, where 
4 cells is the length of the columella and QC. By cutting off 



the root tip of k cells in length, PM length in the decapitated 
root {PM') is N - k cells. 

Variation in k results in different auxin and PIN expression 
dynamics (Fig. 3E). If k was less than 12 cells (PM' > 80 %, 
which means that the decapitated root retained about 80 % of 
the intact PM) the stationary solution had an auxin 
maximum at a certain distance from the root end. The stable 
auxin maximum formed under the reflected flow mechanism 
as described in Mironova et al. (2010). In particular, the 
changes in auxin concentrations calculated in the model fit 
well to changes in DR5 activity detected in vivo during 
RAM regeneration (Fig. 3D). If k was varied between 12 
cells and 20 cells (50 % < PM' < 80 %) damped osciUations 
of intracellular auxin concentrations with several auxin 
maxima appeared in the model solution with a high- amplitude 
zone in the middle of the proximal meristem (Supplementary 
Data Video 2). Finally, if k was greater than 20 cells 
(PM'< 50%; more than half of the intact PM was 
removed), a stationary auxin distribution with the maximum 
in the most shootward cell of the PM (N') was formed, 
decreasing gradually toward the cut edge. These numerical 
results agree well with the experimental data (Fig. 3E) and 
thus explain most of the traits in the RAM regeneration 
processes (see Discussion). 

In silico experiments on excised root tips. To simulate auxin 
distribution in the excised root tip we switched off the shoot- 
derived acropetal auxin flow in the model. In this case, the 




Fig. 4. In vivo and in silico experiments on the root tip excision. (A-D) Time series of the DR5::GFP-expressing excised root tip. (E) Dynamics of DR5 level 
along the central axis of the individual excised root from A-D. (F-I) Model solutions without biosynthesis obtained at different time points as the model 
approached the stationary state. (J) Auxin distribution in the central cylinder layer j = 4 of the solutions (F-I). (K-N) Simulation results of the model with 
auxin biosynthesis fitted to the experimentally observed data (A-D). (O) Auxin distribution in layer 7 = 4 of the model with auxin biosynthesis. 
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Fig . 5 . In silico analysis of maintenance of the auxin maximum in the excised root tips under different conditions. (A) Time of maintenance of the auxin maximum 
in the model with or without reflux. (B) Stability of the auxin maximum under different rates of auxin degradation. (C) Estimated rates of wound-induced auxin 
synthesis for excised root tips of different lengths. Solid line: model prediction. Dashed black line: sectional area of the wound, ^ = 10. It is believed that the 
actual rates of wound-induced auxin synthesis must correlate positively with the area of the wound. 



model calculation was started from the stationary solution with 
the auxin maximum in the fourth cells (Fig. 2D) with a = 0 
(other parameters remained unchanged; see Supplementary 
Data Table S2). We also simulated auxin distribution in the 
excised tips of different lengths. A set of models was created 
with cellular layouts that consist of a similar number of rows 
M = 8, but a different number of cells in a row N = 10, 20 
and 35 (Supplementary Data Table S3). In the model solutions 
we observed a gradual decrease of auxin concentration in the 
root end (Fig. 4F-J), consistent with the experimentally 
observed decrease in DR5 activity over the first 2 d 
(Fig. 4A-C). However, in silico these dynamics were followed 
by the disappearance of the auxin maximum (Fig. 4I-J; 
Supplementary Data Video S3). This disappearance was only 
rarely observed in vivo during the first 5 d after excision (see 
following subsection). Despite this disagreement with the 
experimental data, due to the reflux in the epidermis the dual- 
mechanism model kept the auxin maximum in the excised root 
tip for a longer time as compared with the model without 
reflux (Fig. 5A), as was demonstrated by varying the constant 
for shootward active transport rate (^o,iat,piN2) (Supplementary 
Data Table S3). 

The role of auxin biosynthesis in maintenance of the auxin 
distribution pattern in the root 

In the QC of the excised tips, DR5 activity was maintained 
for a long period (even increasing in the isolated radicles) 
despite the absence of a shoot-derived auxin supply. An intri- 
guing example of the persistent auxin response was observed 
in the shortest excised tips, which consist only of promeristem 
and the root cap (Supplementary Data Fig. SI). Here DR5 
activity in the tip is maintained for a long time despite the 
absence of both the reverse fountain and the reflected flow 
mechanism. Hydrolysis of auxin conjugates or auxin synthesis 
(wounding-induced or ubiquitous) in the excised root tips may 
cause the prolonged maintenance of the auxin maximum. We 
thus formed three hypotheses for the mechanisms that may 
prolong maintenance of the auxin maximum, tested them in 
silico and chose the most preferable. 

Auxin conjugation is one of the pathways for auxin inactiva- 
tion (reviewed by Ljung et al, 2002; Woodward and Bartel, 
2005). By decreasing the value for the auxin degradation 



rate /^d,a we can simulate hydrolysis of auxin conjugates in 
the excised root tip. The decrease in K^^^ prolonged mainten- 
ance of the auxin maximum in the tip (Fig. 5B). However, 
auxin conjugates are the finite resource, which will eventually 
run out without auxin biosynthesis. Thus, we conclude that 
hydrolysis of auxin conjugate may explain the short-lived 
maintenance of auxin distribution but cannot explain the 
prolonged DR5 activity in the excised root tip. 

Auxin synthesis can be rapidly induced by wounding, dir- 
ectly at wounded sites (Sztein et al, 2002; LeClere et al, 
2010). To test whether wounding-induced auxin synthesis is 
able to maintain the auxin maximum in the excised root tips, 
a number of modified models were created wherein auxin syn- 
thesis (eqn 2b) occurred in the last {N\h) cells of the root tip 
for N= 10, 20 and 35 (see Supplementary Data). The rates 
of wound-induced auxin synthesis (/Tg a) cannot be measured 
directly so we estimated them from the stationary model solu- 
tions. The rates were estimated for excised roots of different 
lengths, by requiring that in each case the auxin maximum is 
maintained in stationary solutions despite the absence of 
shoot-derived auxin flow (Fig. 5C, solid curve). We can 
assume that the rates of wound-induced synthesis correlate 
with the cut sectional area, which in turn increases with the 
distance from the root end (Fig. 5C, dashed curve). 
However, the length of the excised root tips did not correlate 
with the rates of wound-induced auxin synthesis estimated in 
the model and estimated on the basis of the cut sectional 
area (Fig 5C; Supplementary Data Table S3). Thus, although 
wound-induced auxin synthesis has an impact on auxin 
supply in the excised tip, this effect is not sufficient to 
produce the experimentally observed data. 

Cells of all types in the root tip have the capacity for auxin 
synthesis and this capacity contributes substantially to robust 
auxin gradients in the root tip (Ljung et al, 2005; Petersson 
et al, 2009). To test the role of the ubiquitous auxin synthesis 
in silico we created a modified model in which all cells ubiqui- 
tously produce auxin at a low rate (eqn 2b). The model 
parameters for auxin synthesis and shoot-derived auxin 
flow to the RAM were adjusted in this model to obtain the 
auxin maximum in the fourth cells of the central cylinder 
(Supplementary Data Table S3). Simulation of the root tip 
cut was then performed as described in the previous section 
for N = 10, 20 and 35 cells. Unlike the case of wound-induced 
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synthesis, in this case, the rates of auxin synthesis were the 
same for root tips of different lengths (Supplementary Data 
Table S3). In the model solutions we observed an overall de- 
crease in auxin content at the root end, although the concentra- 
tion maximum was still located in the 4th cells (Fig. 4K-0; 
Supplementary Data Video S4). Moreover, the model with ubi- 
quitous auxin synthesis equilibrates faster while simulating 
auxin distribution in intact and decapitated root, a possible 
adaptive feature of the system (data not shown). 

Our simulations demonstrate that compared with hydrolysis 
of auxin conjugates and wound-induced synthesis, ubiquitous 
auxin biosynthesis is the preferable mechanism for mainten- 
ance of the auxin maximum in the excised root tips. Thus, 
auxin biosynthesis appeared to be a third force of equal im- 
portance to the reverse fountain and reflected flow mechanisms 
for auxin patterning in the RAM. 

DISCUSSION 

Experiments on cut roots to study meristem regeneration and 
maintenance are classical in plant developmental biology 
(Torrey, 1957; Feldman, 1980; Ivanov, 2007). The molecular 
mechanisms governing RAM development and maintenance 
are orchestrated by the hormone auxin (Sabatini et al, 1999; 
Blilou et al, 2005; Tanaka et al, 2006; Peer et al, 2011). 
Auxin can be regarded as a morphogen, which is involved in 
establishment of patterns of cell identities within developing 
plant organs (Benkova et al, 2009; Wolpert, 2011). A 
unique feature of auxin is its tightly regulated cell-to-cell 
polar transport, which allows auxin to generate gradients and 
maxima in a particular tissue and contribute to tissue pattern- 
ing in a dose-dependent manner (Peer et al, 2011). 
Redistribution of the auxin gradient provides positional infor- 
mation for root regeneration (Xu et al, 2006; Sena et al, 2009; 
Mironova et al, 2010). Auxin regulates expression of its own 
transporters. At relatively low level, auxin upregulates PIN 
family transporters (Vieten et al, 2005). Polar localization 
of PIN proteins on the cell membrane may also be upregulated 
by auxin (Paciorek et al, 2005; Sauer et al, 2006). A negative 
effect of high auxin doses to PIN expression and polarization 
has also been shown by Vieten et al (2005). 

An understanding of the complex interplay of auxin trans- 
portation mechanisms requires the application of sophisticated 
methods. Auxin distribution has been simulated in mathemat- 
ical models predicting significant mechanisms for phyllotaxis 
(Jonsson et al, 2006; Sahlin et al, 2009), tissue polarization 
(Wabnik et al, 2010), leaf venation (Feugier et al, 2006; 
Bayer et al, 2009), shoot branching (Prusinkiewicz et al, 
2009) and root development (Swarup et al, 2005; 
Grieneisen et al, 2007; Laskowski et al, 2008; Mironova 
et al, 2010; Band and King, 2012). Two conceptual models 
have been proposed to explain the mechanism of auxin 
maximum formation in the RAM: the reverse fountain and 
the reflected flow. The reverse fountain (Grieneisen et al, 
2007) is a structural mechanism which implies that the RAM 
structure determines the auxin distribution pattern. The 
reflected flow (Mironova et al, 2010) is a morphogenetic 
mechanism which establishes positional information for 
RAM patterning. It seems intuitively likely that both structural 
and morphogenetic mechanisms act in the root, providing for 



self-organization and subsequent self-maintenance of the 
RAM. Here we combine the two mechanisms in one model 
to test whether they cooperate and what processes in RAM 
development they jointly explain. 

The mechanisms of the reverse fountain and reflected flow 
naturally cooperate in intact roots 

For the dual-mechanism model we used the same cellular 
layout as in the reflected flow model (Fig. lA; Mironova 
et al, 2010). Only two cell types were considered: central cy- 
linder cells and epidermal cells. PINl was the auxin transport- 
er in the reflected flow model. In the dual-mechanism model 
two additional auxin transporters, PIN2 and PIN3, were 
included to mediate shootward and lateral (outward/inward) 
auxin transport. The dual-mechanism model integrates 
up-to-date experimental data on auxin-regulated PIN1-PIN3 
expression (Fig. IB-D). The model parameters for PIN ex- 
pression were estimated so that the effective auxin concentra- 
tions for PIN regulation calculated in silico match those 
reported in vivo by Vieten et al (2005; see Supplementary 
Data). In the model solution we observed formation of PINl, 
PIN2 and PIN3 expression domains together with auxin distri- 
bution pattern, both agreeing well with the experimental data 
(Fig. 2). 

We suggest that the two mechanisms may naturally 
cooperate in vivo due to auxin-regulated expression of its 
own transporters. By using model analysis we estimated the 
concentration of auxin, PINl, PIN2 and PIN3 in different 
RAM tissues. The estimated auxin concentrations (0-4-1 |jlm 
in a vascular cell; 10-11 |jlm in the QC and 8- 10 |jlm in a colu- 
mella cell) agree well with the following constraints from in- 
dependent data. DR5 reporter expression pattern was 
demonstrated to reflect free-auxin concentrations between 
10"^ and 10""^ m (Ulmasov et al, 1997; Sabatini et al, 
1999). Also, Petersson et al (2009) showed that the concentra- 
tion of intracellular indole-3 -acetic acid in the RAM is in the 
range 1-50 |jlm. 

In numerical experiments simulating self-organization of 
auxin distribution in an initially unspecialized tissue, we 
observed generation of an auxin distribution pattern that 
matches the experimentally known expression of DR5. The 
auxin maximum peak in the root tip is formed under the 
reflected flow mechanism, which may operate from the earliest 
stages of root development. Thus, cells of the developing root 
can acquire positional information sufficient for RAM pattern- 
ing. Moreover, in the model solution all necessary features for 
the functioning of the reverse fountain were organized by the 
reflected flow mechanism. So, the further maintenance of the 
auxin maximum at later developmental stages can be provided 
mainly by the reverse fountain structural mechanism. The 
auxin gradient in the PM is also organized by the reverse foun- 
tain. However, the reflected flow mechanism does not dis- 
appear even after formation of the RAM structure. The 
morphogenetic mechanism continues to function, revealing 
itself if the RAM structure is disrupted or the environment 
changes. For example, the mechanism is crucial for RAM 
re- specification after NPA (naphthylphthalamic acid) treatment 
(Mironova et al, 2010) or regeneration of the RAM in 
decapitated root. 
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The next step for model testing was a validation experiment 
to test whether the model would fit the evidence and provide 
reliable predictions. For the dual-mechanism model, experi- 
ments on root cut are germane. Performed from the beginning 
of the 19th century in many labs, experiments on the regener- 
ation of the decapitated root and on the viability of the excised 
root tips outlined phenomenologically the steps occurring in 
the two root fragments. The advance in this class of experi- 
ment here was in using DR5::GFP A. thaliana plants that 
reveal an auxin response gradient in addition to classical 
phenomenology. 

The mechanism for auxin distribution self-organization 
in the proximal meristem of decapitated roots 

Previously, two main types of root restoration in decapitated 
roots were reported. If the removed part is sufficiently large 
(proximal cuts, for example 1 mm for Vicia faba or 700- 
750 |jLm for Zea mays) then lateral roots formed close to the 
cut edge (Feldman, 1980). However, removal of a smaller 
length of root tip (distal cuts, for example 450 ixm for 
Z. mays) led to regeneration of the root tip, aligned with the 
main axis of the root. Distal cuts were carefully investigated 
by Sena et al (2009) and it was shown that the region of 
high regeneration competence in the A. thaliana root was 
located in the distal PM (up to 130 ixm from the root end), 
the frequency of regeneration decreased by a factor of ten 
more proximally (at 200 ixm from the root end) and dropped 
to zero at 270 ixm. Feldman (1980) showed that in distal 
cuts the part of the PM left in the decapitated root should 
exceed some critical size for root tip regeneration to occur. 
It was also shown that regeneration of the apex arises from 
the actively proliferating PM cells. 

In our in silico and in vivo experiments we found three types 
of response (Fig. 3E). No regeneration occurred in vivo if the 
decapitated root had less than half of the initial PM. 
Decapitated roots containing approx. 50-80% of the PM 
first formed a callus-like outgrowth from which a new root 
tip appeared. The meristems of decapitated roots containing 
more than 80 % of the PM regenerated similarly to observa- 
tions for QC laser ablation (Xu et al, 2006). 

Based on our numerical experiments we propose a mechan- 
ism explaining the existence of these three different scenarios 
of new root tip development and their dependence on the size 
of the PM that was left in the decapitated root. The key role is 
played by shoot-derived auxin flow to the RAM, which was 
not changed by decapitation. If the major part (more than 
half) of the PM was lost, auxin flow from the elongation 
zone was too high to be assimilated in the remaining part of 
the PM so that the cells of PM accumulated auxin and this 
did not facilitate meristem restoration. In contrast, if the 
major part of the PM was left in the decapitated root, the 
auxin transport system (namely, the reflected flow mechan- 
isms) was able to cope with the rootward auxin flow, generat- 
ing an auxin maximum at a distance from the root edge. In 
intermediate cases, we observed in silico damped auxin oscil- 
lations with several unstable maxima caused by disturbance of 
the auxin traffic. We suggest that the damped oscillations pre- 
determine callus-like outgrowth. In fact, discordant positional 
information due to temporally varying auxin maxima may 



cause outgrowth cell specialization, making these cells differ- 
ent in size and DR5 level. 

Excised root tip as an autonomous organ 

Prolonged maintenance of RAM viability in excised root 
tips has been known for almost a century (reviewed by 
Torrey, 1957). Study of the role of auxin in RAM stem cell 
niche maintenance has become possible in recent years using 
auxin sensors. Grieneisen et al (2007) showed maintenance 
of the auxin maximum in dissected roots from shoots. 

Here we describe for the first time the dynamics of DR5 
activity in excised root tips of different lengths. Taking into 
account that DR5 is a promeristem marker, the experimental 
results provide evidence for prolonged maintenance of the 
stem cell niche in the excised root tips. The auxin response 
maximum is maintained even in the shortest tips lacking any 
PM. 

Grieneisen et al (2007) suggested that the reverse fountain 
mechanism provides for a stable auxin maximum in dissected 
roots. Consequently, we tested whether this feature is applic- 
able to the dual-mechanism model. The model showed that 
the reverse fountain mechanism has a significant impact on 
auxin gradient maintenance in the excised root, which 
depends on the length of the root tip (Fig. 5A). However, in 
numerical simulation the gradual drop in auxin concentration 
causes the eventual disappearance of the auxin maximum. 
This simulation result suggests that the reverse fountain and 
the reflected flow together were not enough to describe the 
whole set of experimental data. Here we explored in silico 
three possible hypotheses for maintenance of the auxin 
maximum in excised root tips: (1) auxin conjugates hydrolysed 
back to the free hormone can act as sources of free auxin; (2) 
wounding-induced auxin synthesis supports the maximum; 
and (3) the existence of a basal auxin synthesis in root cells. 

Our simulations showed that three mechanisms have differ- 
ent impacts on the maintenance of the spatial distribution of 
auxin. As auxin conjugates are a finite resource, their hydroly- 
sis can account only for the maintenance of the auxin pattern 
for a short period. Soeno et al (2010) demonstrated that after 
treatment of seedlings with an inhibitor of auxin biosynthesis, 
endogenous auxin levels in roots decreased to less than 10 % 
of the levels in the controls, indicating that auxin synthesized 
de novo is the main auxin source for roots. 

The root tip has the ability to synthesize auxin (van 
Overbeek, 1939; Feldman, 1980; Ljung et al, 2005; 
Stepanova et al, 2008; Petersson et al, 2009). van Overbeek 
(1939) proposed that the excised root tip is maintained alive 
due to continuing auxin biosynthesis. The question is 
whether it is a wound-induced synthesis or is there a basal cap- 
acity for auxin synthesis in root cells. Our numerical experi- 
ments confirm that wound-induced auxin synthesis is not the 
main contribution for auxin supply to the excised root tip. 
For this to occur, the rates of auxin synthesis at the wounded 
sites would have to be higher for small wounds and lower 
for large wounds. 

Auxin synthesis was found in both meristematic and non- 
meristematic parts of the root, indicating that there is basal 
auxin synthesis in root cells (Ljung et al, 2005). Our simula- 
tions showed that this basal auxin synthesis together with the 
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reverse fountain allows the auxin gradient to form, with a 
maximum at the 4th -5th cell from the root end in the 
excised root tips. The more PM was left in the excised root 
tip, the more impact the reverse fountain has in maintaining 
the auxin gradient. Basal auxin synthesis can explain mainten- 
ance of the auxin maximum in very short excised roots, which 
contain only QC and the root cap (neither the reflected flow 
nor the reverse fountain mechanism work in excised root 
tips of this type). In our simulations we showed that the 
auxin distribution equilibrates faster if we take into account 
auxin biosynthesis. This can be an adaptive feature for RAM 
development. Rowntree and Morris (1979) showed that ^^C 
applied to apical pea buds accumulated in the lateral root prim- 
ordia but not in the root tip of the primary root nor, in older 
plants, in the tips of the lateral roots. These data support the 
idea that the root tips are autonomous from shoot-derived root- 
ward auxin flow because of de novo auxin synthesis therein. 
This explains why the auxin gradient persists in the root tip 
of the excised root. 

Several recent studies have revealed auxin biosynthesis 
mechanisms in root (Ljung et al, 2005; Stepanova et al, 
2008). Here we suggest that auxin biosynthesis is a third 
main force in auxin gradient formation along with the reflected 
flow and reverse fountain transportation mechanisms. We 
showed that the mechanisms cooperate, providing for a 
semi-autonomous RAM with the ability to self-regenerate. 

SUPPLEMENTARY DATA 

Supplementary data are available online at www. aob. oxford- 
journals. org and include: a description of the 2-D dual mech- 
anism model; the dual-mechanism model reduction to the 
reverse fountain model; parameter estimation (incorporating 
Table SI: parameter estimation for the dual-mechanism 
model, and Table S2: the basic set of model parameters for 
the 8 X 20 2-D dual-mechanism model); simulation of auxin 
distribution in the excised root tips (incorporating Table S3: 
parameter values for the model's modifications); and 
Fig. SI: auxin response maintenance in the shortest excised 
root tips. Video SI: simulation of auxin patterning in the dual- 
mechanism model. PINl, PIN2 and PIN3 expression domains 
and the pattern of auxin distribution are self-organized in the 
in silico root. Video S2: oscillations in auxin concentration 
simulated in the decapitated root may predetermine callus-like 
outgrowth. The model was run from the reduced cellular layout 
with the same intensity of auxin influx to the RAM. Video S3: 
changes in auxin distribution in model simulation without 
auxin biosynthesis. Simulation in the root tips of different 
lengths is presented. Video S4: comparison of the maintenance 
of the auxin maximum in excised root tips with and without 
auxin biosynthesis. 
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